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Abstract 


Supersymmetric 5'O(10) grand unified models with renormalizable 
Yukawa couplings involving only 10 and 126 Higgs fields have been 
shown to realize the fermion masses and mixings economically. In pre¬ 
vious works, the sum rule of the fermion mass matrices are given by 
inputting the quark matrices, and the neutrino mixings are predicted in 
this framework. Now the three neutrino mixings have been measured, 
and in this paper, we give the sum rule by inputting the lepton mass 
matrices, which makes clear certain features of the solution, especially 
if the vacuum expectation values of 126 -|- 126 {vr) are large and the 
right-handed neutrinos are heavy. We perform the analyses to fit the 
fermion masses and mixings using the sum rule. In previous works, the 
best fit appears at vr ~ 10^^ GeV, and the fit at the large vr scale (~ 10^® 
GeV) has been less investigated. Our expression of the sum rule has a 
benefit to understand the flavor structure in the large vr solution. Using 
the fit results, we perform the calculation of the /i —)■ ey process and the 
electric dipole moment of electron, and the importance of vr dependence 
emerges in low energy phenomena. We also show the prediction of the 
GP phase in the neutrino oscillations, which can be tested in the near 
future. 


PACS. 12.10.-g, 12.10.Dm, 12.10.Kt 




1 Introduction 


The Higgs discovery at the LHC [T] opens a new era to nnderstand the fermion masses and 
the mixings, which are generated by the Ynkawa interaction to the Higgs boson. Indeed, most 
of the parameters in the standard model (SM) he in the Ynkawa conpling matrices. As the 
remaining parameters in the flavor sector in the SM, the nentrino 13-mixing angle has been 
measnred by nsing the short baseline nentrino oscillations from the reactor nentrino [2], and 
the CP phase in the three-flavor nentrino oscillations (the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) phase) is expected to be measnred accnrately in the long baseline oscillations [3]. 
The preferable region has been already obtained by the global data analysis for the nentrino 
oscillations HEIE]. Do the accnrate measnrements of the parameters in the SM mean the 
end of the story? The answer is of conrse No! Even after the parameters are accnrately 
measnred, they do not £x the strnctnre of the Ynkawa matrices becanse of the existence of 
the nnphysical flavor rotation. The particle physicists claim that there mnst be new physics 
beyond the SM from varions points of view. The nnphysical rotation in the SM can be physical 
in models beyond the SM. In the era that the SM parameters in the flavor sector are accnrately 
measnred, we shonld, thns, stndy the strnctnre of the Ynkawa matrices which can inflnence 
observables beyond the SM to test the varions flavor models. 

The difference of the generation mixings between qnark and lepton sectors is one of the 
major issnes to nnderstand the strnctnre in the Ynkawa conpling matrices in the nnihed flavor 
pictnre. The S'O(IO) grand nnihed theory (GUT) provides a promising framework to nnify 
the qnarks and leptons, becanse the entire SM matter contents of each generation (inclnding 
a right-handed nentrino) can be nnihed in a single irredncible representation, 16 . Particnlar 
attention has been paid to the renormalizable minimal AO (10) model, where two Higgs mnlti- 
plets { 10 © 126 } are ntilized for the Ynkawa conplings with the matter representation [7]. The 
conplings to the 10 and 126 Higgs helds can reprodnce realistic charged fermion mass matrices 
nsing their phases thoronghly [8]. Qnalitatively, the smallness of the generation mixings in the 
qnark sector can be explained by the left-right symmetry which is a snbgronp of the AO(IO) 
symmetry, while the nentrino mixings are not necessarily small and large mixing angles for 
the solar and atmospheric nentrino oscillations can be obtained natnrally in this framework. 
Since the fermion mass matrices are given by the linear combination of two symmetric matri¬ 
ces (mnltiplied by donblet Higgs mixings), one can obtain the snm rnle of the fermion mass 
matrices. Becanse the nnmber of the physical parameters in the minimal model is restricted, 
algebraic predictions of the SM parameters can be obtained nsing the snm rnle. Actnally, the 
qnantitative prediction of nentrino oscillation parameters has been discnssed by inpntting the 
qnark masses and mixings [a [10]. The model predicts that the nentrino 13-mixing angle is 
non-zero, and the predicted valne is consistent with the measnred valne. We now repeat the 
qnestion in the previons paragraph. When the mixings and the CP phases in the nentrino 


1 




oscillations are measured accurately, does it mean the end of the study of the minimal SO(10) 
just by checking with the predictions after the number-crunching hts? 

The supersymmetric (SUSY) theory is one of the most attractive candidates to construct 
a model beyond the SM, and can explain the large hierarchy between the weak scale and the 
GUT scale. In SUSY models, the structure of the Yukawa matrices can influence the low 
energy phenomena if the SUSY particle masses are less than around a few TeV. Even after the 
SM parameters are htted (within the experimental errors), there can be physical degrees of 
freedom to influence the observables, though they also depends on the SUSY particle spectrum. 
Therefore, it is important not only to £t the SM parameters, but also to understand the flavor 
structure consistent with the experimental data in the format applicable to extended models. 
From such a point of view, in this paper, we revisit the hts of the fermion masses and mixings 
in the minimal SUSY S'O(IO) model. We hrst describe the sum rule of the fermion mass 
matrices in terms of the charged lepton mass matrix M^, and seesaw neutrino mass matrix 
M.V as inputs, and express the quark mass matrices by them. In fact, the quark masses and 
mixings can receive large radiative corrections at low energy threshold in the supersymmetric 
version of the model, and there can be large ambiguities. Whereas, the masses and mixings 
in the lepton sector are now more accurately measured and less ambiguous. Therefore, it is 
essential to investigate the properties of the solution inputting the parameters in the lepton 
sector. Using the expression by inputting the parameters in lepton sector, we perform a 
analysis to fit the quark masses and mixings. Surely, the result must be the same whichever 
input is chosen in the expression. Our claim is that the algebraic expression can be very helpful 
to understand the structure of the Yukawa matrices which reproduces the experimental data. 

The scale of 126 -I- 126 vacuum expectation value (vev), vr, is important for breaking 
the rank-5 S'O (10) symmetry down to the rank-4 SM gauge symmetry. The scale is also 
important to know how the light neutrino mass is generated by the seesaw mechanism m- 
In the S'O(IO) model with 126 Higgs coupling, not only the right-handed neutrino Majorana 
mass for type I seesaw, but also the SU{2)r triplet contribution to the light neutrino mass is 
generated known as type II seesaw na. In GUT models, an intermediate scale often appears, 
and it can implement the proper size of the light neutrino masses. In fact, the best fit in the 
minimal model is obtained when vr ~ 10^^ GeV, as was reported in Refs. [T3| [14]. However, 
if Vr is about 10^^ GeV in SUSY models, the decomposed representations in 126 -|- 126 he 
around the scale, and they are harmful to the gauge coupling evolution since gauge couplings 
may blow up before the unification. Being motivated by this fact, we feature the solutions for 
Vr being around 10^® GeV. In order to construct the solutions in which the type I seesaw part 
contributes the neutrino masses, the Yukawa coupling matrix / for 126 Higgs coupling should 
be nearly singular and the elements of f~^ are enlarged. Such a situation is not obvious if the 
quark masses and mixings are inputs. Gontrary to the situation, in our expression in terms of 
the parameters in the lepton sector, it is easy to decode the structure of the singular matrix. 
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We show the results of the fit by changing vr, and also show the result of type I seesaw, 
where SU (2) l triplet contribution is negligible. From the ht results, we point out that future 
accurate measurements of the neutrino 23-mixing angle and the PMNS phase play a decisive 
role to determine the Yukawa structure in the solutions for vr > 10^® GeV. This feature can 
be algebraically explained in our expression of the sum rule. 

In the minimal S'O(IO) model, the Dirac neutrino coupling and the Majorana neutrino 
mass matrix are predictable depending on the scale vr. As a consequence, by investigating 
the lepton flavor violation (LFV), it is possible to hnd a footprint to know which scenario 
of the neutrino mass generation is chosen. We will perform the calculation of the branching 
fraction of /i ^ ey decay and the electric dipole moment (EDM) of the electron, using the 
prediction from the analysis in the model and discuss the vr dependence of the predictions. 

This paper is organized as follows: In section 2, we introduce the sum rule of the fermion 
mass matrices in the 5'O(10) model with the minimal Yukawa interaction. In section 3, the 
expression of quark mass matrices in terms of the mass and mixing parameters in the lepton 
sector is derived. We examine the properties of the solutions expressed by the parameters in 
the lepton sector, and study the vr dependence of the solutions. We perform the numerical 
analyses to £t the fermion masses and mixings in section 4. The importance of the predictions 
of the phase in neutrino oscillations is illustrated. In section 5, we apply our £t results to 
the calculations of branching ratio of /i —)■ ey process. We also mention the vr dependence of 
the results and study the importance to calculate the electron EDM. We discuss the possible 
modihcation of the Yukawa couplings by the SUSY threshold corrections. Section 6 is devoted 
to conclusions and discussions. In the appendix, we express the square root matrix, which is 
necessary to express the solution in terms of the parameters of the lepton sector. 

2 Sum rule of the fermion mass matrices in the minimal 
*5'O(10) model 

The SUSY GUT can provide the most promising framework to incorporate the vast data 
systematically and consistently. Among many candidates, S'O(IO) [l5] is the smallest simple 
gauge group under which the entire SM matter contents of each generation are unihed into a 
single anomaly-free irreducible representation, 16. The 16-dimensional spinor representation 
in S'O(IO) includes the right-handed neutrino and no other exotic matter particles. 

The minimal S'O(IO) models is dehned as that in which the Yukawa interaction is minimal, 

^ In this paper, we define the minimality of the S'O(IO) model on the Yukawa interaction: The matter 
representations couple to only 10 and 126 Higgs representations, and the number of parameters in the Yukawa 
couplings is minimal (without a symmetry). People often impose the minimality of the model even on the Higgs 
contents, 10 -|- 126 + 126 -I- 210, in which the number of parameters in the Higgs interactions is minimal (in 
addition to the minimality in the Yukawa interactions). In this case, however, the SUSY version of the model 
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namely, only H : 10 and A : 126 Higgs representations couple to the fermions "0 : 16 by 
renormalizable interaction: 

Wy = ( 2 . 1 ) 

Due to the S'O(IO) algebra, the coupling matrices are symmetric, hj^ = hji and = fj/. 

The Higgs superpotential is investigated to see the detail pattern of symmetry breaking 
from S'O(IO) down to the SM gauge group [HI [19]. In these analyses, the criterion of the 
renormalizability plays an essential role not only to reduce the number of parameters but also 
to construct a model without the vevs of H — L = ±1 direction by IQh + 16^ representations 
to break S'O(IO) symmetry. Actually, the vev oi B — L = ±2 direction by 126 + 126 helds 
reduces the rank of the gauge symmetry in the renormalizable model, in which the i?-parity is 
automatically conserved in the minimal SUSY standard model (MSSM) vacua. The extension 
to include 120 Higgs representation without breaking the renormalizability has been also con¬ 
sidered [20]. The 126 Higgs representation includes both SU{2)l and S'f/(2)j^ triplets, and the 
Yukawa interaction to 126 can generate both left- and right-handed Majorana neutrino mass 
matrices. Therefore, in the framework of the minimal S'O(IO), depending on the symmetry 
breaking vacua, the light neutrino mass matrix can be obtained [TOl IT3] in both type I and H 
seesaw mechanism mn 

The Yukawa coupling (after S'O(IO) symmetry is broken down to the SM) is given as 
follows: 


Yu = h + r2f, 

Yd = ri{h + f), 

Ye = nih-3f), 

Y, = h- 3r2f, (2.2) 


where ri and r 2 depend on the Higgs mixing (doublet Higgs mixing in 10 and 126 ), and h 
and / are original Yukawa matrices h and f multiplied by Higgs mixing^: 


h = Unh, 


/ 


Ul2 . 
x/3n ’ 




n = 


11 


Ui 


Ui 


r2 = ri 


13 


11 


U 


12 


(2.4) 


suffers some delicate mismatches between all the fermion data fittings appeared in the recently developed 
measurements and the self-consistency of the gauge coupling unification [14) . and the framework suggests a 
non-SUSY model (or split SUSY) [THldT]. Since we focus on the prediction on the Yukawa structure which 
can communicate with the low energy phenomena, we impose the minimality just on the Yukawa interactions. 
See section 6 for discussion. 

^ The unitary matrices U and V are the diagonalizing matrices of the doublet Higgs matrix Mdoubiet: 
UMdoubietY^ is diagonal, 

— ^doublet = {Hd°’ Ad, Ad, <l>d)Mdoublet(H^°, A„, A„, <i)„)’^. (2.3) 

The lightest linear combinations of the doublets will be the MSSM Higgs doublets. The detail can be found 
in Ref. [21] (and [19] in different conventions). 
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The charged fermion masses are obtained as 


Mu = YuVu, Md = YdVd, Me = YeVd, 


M° = ni.. 


( 2 . 5 ) 


where Vu and Vd are the VEVs of up- and down-type Higgs helds. We obtain the relation of 
the mass matrices as 


Md = Me + -^F = rMu + F, (2.6) 

1 - r2 

rMf = Me + 3F, (2.7) 


where 


- .o 

r = ri — = ri cot p, 
Vu 


and the matrix F, which is proportional to the 126 Higgs coupling matrix, is 


F = r(l - r 2 )fvu 


( 2 . 8 ) 


(2.9) 


Roughly, we obtain r ~ nih/mt. Surely, these mass relations are realized at GUT scale, and 
the evolution via renormalization group equations (RGBs) is considered to ht them with the 
low energy data of quark-lepton phenomena including the neutrino oscillations [9]. 

The right-handed Majorana neutrino mass matrix is obtained as 


Mr = V2ivR, (2.10) 

where vr is a VEV of 126 . Practically, we denote 

Mr = CrVr/, ( 2 . 11 ) 

where 

cr = V6^ = x/6^, (2.12) 

Ui2 ri3 

for the current notation. Because U 12 and U 13 are components of the diagonalization unitary 
matrix for doublet Higgs helds, cr has a minimal value. The size of cr is related to the size 
of original f coupling, which will be important to derive the GUT scale threshold correction 
for havor violation. One can rewrite as 


Mr = 


CrVr 


-F. 


r(l - r 2 )vu 

The seesaw neutrino mass matrix can be written as mm 

M. = ML-M^M],\M^f, 


(2.13) 


(2.14) 
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where Ml is the left-handed neutrino Majorana mass which comes from SU{2)l triplet cou¬ 
pling |22], UAl- In S'O(IO) model, the 126 Higgs also includes the SU{2)l triplet and 
the Yukawa coupling generates both Ml and YTr. Therefore, Ml is also proportional to the 
coupling matrix / and we denote 

Ml = clVlJ. (2.15) 

In the S'O(IO) model, if there is 210 or 54 Higgs represent at iorE, the VEV of Al, {Al) = vl, 
can be obtained as n^gak/^A, where M^ is the mass of the SU{2)l triplet. 

For convenience, we reparameterize Mli and Ml as 

Mr = -^F X 10'°, (2.16) 

Ml = 6Fx 10"'°, (2.17) 

where 

R= X 10'°, (5 = ——X 10'°. (2.18) 

r crVr r(l - r 2 )Vu 

The factor 10'° is attached for a numerical convenience. We will treat i? as a parameter to 
specify the scale of vr implicitly. For example, if vr ~ 10'° GeV, one obtains R 0(1) 
roughly. For vr ~ 10'° GeV, one hnds R ~ 0(10“°). 

We remark on the naive statement on the scale vr to £t the fermion masses and mixings, 
which can be found in the previous works. One may think that the 0.05 eV neutrino mass can 
be easily obtained even if vr ~ 10'° GeV because the scale of Mr = crVr/ can be naturally 
0(10'“') GeV for the component of / to be 0(10“^). However, such a naive thought is true 
only if the atmospheric mixing is small. In fact, in order to reproduce a large atmospheric 
mixing, not only {A4u)33 but also (Aiu )23 should be 0(0.01) eV. In the Mg-diagonal basis, 
one obtains (MeF“'Me )23 (MeF“'Me) 33 , and thus, to obtain a large 23-mixing, one needs 

a cancellation in (3,3) element of and the naive thought above does not work. To obtain 
the large atmospheric mixing, {Mef~^Me) 23 /{cRVR) has to be 0(0.01) eV. Because the size 
of / is 0(0.01) due to the htting of quark masses and mixings, the naive size of the vr 
is O(10'°) — 0(10'^) GeV, which corresponds to i? ~ 0(1)- Oiie can find a solution for 
Vr ~ O(10'°) GeV (i.e. R ~ 10“°), only if the matrix / is nearly singular and /“' is enhanced 
(compared to its naive size from each component). 

We note that the solution for the nearly singular / matrix is also found in Ref. [H] (they call 
the solution as “Mixed'”). However, if one fits the charged fermion masses, there is no reason 
that / is close to a singular matrix, and it is not easy to study the property of the solution. 
We describe the neutrino mass and mixing parameters to be inputs, and construct a solution 
where the situation of small det / can be seen explicitly. The solution of nearly singular Mr 

^ When 54 is employed to generate the triplet contribution to the neutrino masses [23], the S'O(IO) breaking 
vacua is modified, and r 2 parameter can be independent to the vacua. 
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matrix is applicable even to the non-minimal models. We believe that the description of nearly 
singular Mr is useful to study how the light neutrino mass scale is obtained in seesaw models. 

3 Property of the solution by inputting the lepton pa¬ 
rameters 

We have described the general setup in the minimal S'O(IO) model. In the previous works, 
the htting of the fermion masses and mixings has been performed by inputting the quark 
masses and mixings, and the neutrino mixings are predicted in the framework. At present, 
the lepton parameters are more accurate rather than the quark ones, and besides, a large 
threshold correction are expected in the quark sector in SUSY models. In that sense, it is 
better to perform the htting by inputting the parameters in lepton sector, and outputting the 
quark parameter. The formulation is presented in not only such a practical purpose, but also 
to make clear the property of the solution with vr ~ 10^® GeV, which is main concern in this 
paper. 

Using the relation, rM^ = + 3F, we obtain 

X 10^° = {6F - R{Me + 3F)F-\Me + 3F)) 

= {{6-9R)F-6RMe-RMeF-^Me) . (3.1) 

This equation can be rewritten as a quadratic equation in terms of and we 

obtain 

F = - (k-QR1 + VK‘^ - 12RK + AdR l) (3.2) 

18R -25 \ ) 

where 

K = M^ = M^x 10^°. (3.3) 

The square root matrix is dehned in the Appendix. Using the expression of F matrix, the 
quark mass matrices, and M^, can be given as a function of charged lepton mass matrix 
Mg and the light neutrino mass matrix ATj/. 

Now, by using the formula given in Appendix, let us express the matrix 

K-6R1 + ^(K-6R1)^ + B1, (3.4) 

where B = 45R — 36R^. Surely because K and (K — 6R1)^ + B1 can be diagonalized 
simultaneously, we can easily obtain 

K-6R1 + ^/{K - QRiy + Bl = '^(^Xi-QR + - 6R)^ + B^ A^, (3.5) 
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where Xi {i = 1,2,3) are eigenvalues of matrix K, Si are signs for the square roots, and A* 
are given in Appendix. By the dehnition of the matrix K, K = —Me M.yMe (where 
My = My X 10^°), we obtain 




1 

(Ai — A2 )(Ai — As) 


(^MyM^ ^My + (A 2 + Xs)My + \2X2,M^ , 


(3.6) 


and similarly for A2_3. 

In the limit of i? —)■ 0 (i.e. vr —)■ 00), we obtain 

F (k + ^/lo) (3.7) 

2o V / 

So, let us evaluate K + \fK^. One can write down 

(1 + si)Ai A 

K + ^ = V\ (1 + S2)A2 y = 5^(1 + s,)A,A,. (3.8) 

(1 + S3)A3 ) i 


For example, suppose Si = —1 and 82,3 = 1 (because the square root of a complex number has 
two branches, we choose a sign convention \/~M = A to dehne Sj), one hnds 

K + \/l0 = 2 ^ A*A* = 2K - 2A1A1. (3.9) 

i=2,3 


As a result, we obtain under the choice of Si and in the limit of i? —)■ 0, 


F 


^ [jCiy + AiMi/^AiM^/^j 

TTT-, N / ,- 7 ^ (XiMyM^ ^Mv + (Ai + X^X'^My + A1A2A3 

0 (Ai — A2 )(Ai — A3) V 



(3.10) 

(3.11) 


In the choice of si = —1, S2 = S3 = 1, it can be written as 


/ 0 

K + \/)^ = 21/ A2 

V As 




(3.12) 


Therefore, rank(iF + VK^) = 2 in this case. As a consequence, one of the right-handed 
neutrinos is very light compared to the other two because the right-handed Majorana mass 
matrix is proportional to F. Needless to say, R is hnite (though it is small), and therefore, 
the lightest right-handed neutrino is not massless. 

For a choice of si = S2 = S3 = 1, one obtains 


K + -/lO = 2K, 


( 3 . 13 ) 




and thus, in the limit of i? —)■ 0, 


( 3 . 14 ) 


F ^ Im^. 

0 

Therefore, this choice corresponds to the case where the right-handed neutrinos are decoupled, 
and the triplet contribution is dominant (namely = Ml)- It is known that the fits of 
fermion masses and mixings are not good in the triplet-dominant case [H]. The difference 
compared to the previous case comes from the contribution from which makes 

the hts better. In our description of the solution, the triplet-dominant case is contained in the 
expression just by a choice of the signs of the square root matrix. In fact, as we execute 
analysis in the following sections, the triplet-dominant case does not appear in the fits by 
minimum, because it gives a large value. Therefore, we call the combination of the triplet 
and type I seesaw contributions just as type II solution, which is intrinsic in the S'O(IO) model. 
In the Ref. [ll] , they name the solutions as “Mixed” and “Mixed'”, but those two solutions are 
continuously connected as a function of the parameter R in our description. 

We note that the dependence of R is small in the type II solution for si = — 1, S2 = ss = 1 
for small R (namely, large vr) as one can find from the i? —?■ 0 limit. In type I seesaw (5 = 0), 
the F matrix in the limit of i? ^ 0 can be written as 

F ^ [k + /P) MlF, (3.15) 

and it surely depends on R. For the solution si = S2 = S3 = 1, one obtains 

F^~M.. (3.16) 

For a smaller value of R, the component of F becomes large and it is expected that there is 
no solution. For a choice of Si = S2 = S3 = —1, one obtain^ 

F ~ RMl^^K-^Ml/'^ = -RMeM-^Me. (3.18) 

For i? —)■ 0, F becomes small in this choice and which means there is up-down symmetry 
{Md (X Mu) in the limit, and there is no solution. Therefore, there may no solution for 
"Vr > 0(10^®) GeV in type I. 

Because of the algebraic equation {K—6R l+\/K'^ — 12RK + 4:dRl)(K—6Rl—\/K'^ — 12RK + ASRl) = 
2R{18R — 2(5)1, one obtains 

F = - - -- 6i? 1 - - 12RK + A6R 1)mF2 

18i? — 26 

= 2 RMIF{K -6R1 + - 12RK + 4 :SR 1 )~^mIF, (3.I7) 
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4 Numerical verification 

In the previous section, we have studied the property of the solution for the case of a large 
vr (/?—)■ 0). Now, we verify the property by performing the numerical htting of the fermion 
masses and mixings. In our formulation, the lepton parameters are inputs and the quark 
parameters are outputs. Of course, in order to perform the numerical ht in the practical 
way, we use the observed quark masses and mixings to £t by analyses. Therefore, the £t 
results have to be the same as the one given in the previous works (up to the detail updated 
experimental data) mathematically, and the formulation itself does not provide a better htting 
for sure. As described, the purpose to use the expression by inputting the lepton parameters 
is to make clear the property of the solution depending on the vr, which is important to 
apply the model to the low energy physics. In fact, the htting in the minimal 5'O(10) model 
needs complicate tuning in the parameters (depending on the parameterization) especially in 
the case of a large vr, and one may feel that a better ht can exist by using the tuning of 
the parameters in quark sector. However, in our expression, the solution provides a simple 
behavior for the large vr case, and we can understand how the ht behaves. 

In this section, we demonstrate the numerical ht of the fermion masses and mixings, and 
we check the consistency to the previous works in the literatures. Compared to era when the 
previous works are done in the literatures, the neutrino oscillation experiments enters to a new 
stage to measure the oscillation parameters: not only the measurement of the 13-mixing angle, 
but also the expectation of the precision measurements of 23-mixing angle and PMNS CP 
phase. We hrst describe our parameterization which is suited in our scheme, and next explain 
the method of our analyses including a technical detail. As a purpose to revisit the analyses, 
we show our ht results making clear the importance of the precise measurement of 23-mixing 
angle. Although the ht of the PMNS phase is constrained mildly, the correlation between the 
PMNS phase and 23-mixing angle can provide an important implication to distinguish models 
and the havor structure in the diherent essential parameter regions. 

4.1 Parametrization 

The relation is summarized as 

rM„ = + 

1 - ra 

Mrf = rMu + F, 

F = -^- -My^ (k -6R1 + - 12RK + 

18R-2S V / 


(4.1) 

(4.2) 

(4.3) 

(4.4) 
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We work on the basis where the charged lepton mass matrix Me is diagonal: 

Me = diag. (mee*“%(4.5) 

The light neutrino mass matrix is parametrized as 

= ^diag. (4.6) 

and the unitary matrix U is the PMNS neutrino matrix, and we use the usual convention 
by particle data group (PDG). In the matrix 17, there are three neutrino mixing angles 
^12,^23,^13 and a CP phase ^pmns in neutrino oscillation. As a convention, the mass pa¬ 
rameters, mi, m 2 , m 3 are given as real and positive values. 

As one can hud easily, one of the phases ccg, ci/x, cir,«!, <^2, CI3 can be made to be zero 
without loss of generality, since the relation is covariant under 

^ M, ^ (4.7) 

We choose ai = 0 as a convention. The relation is also covariant under the rephasing: 

R P°‘R, 5 —?> Mu^d,e e (4-8) 

and thus we can choose R to be real (and positive) without loss of generality. We can also 
dehne the parameter r to be real. Therefore, the parameters in this parametrization are the 
following: 


me,m^,mr,mi,m2,m3, 

^65 ^ 2 ? *^35 

6^12, ^23, ^13, ^PMNS, (4.9) 

r2 (complex), r, 

(5(complex), R, 

and there are 21 degrees of freedom in total. 

The input parameters in the lepton sector are 

me, m^, mr, Q 12 , 6 * 13 , 6 * 23 , (5pmns, (4.10) 


and the PMNS phase (5pmns will be treated as output in our analysis. The parameters in the 
quark sector are 

mu, me, mt, md, ms, mb, Vus, Vcb, Vub, (5km, (4.11) 


which we will £t by using the freedom in the model, 
observablefl 


Totally, there are 19 degrees for the 


® Depending on the size of absolute neutrino mass, the neutrino Majorana phases 02 and 03 and the 
neutrino mass mi can be the observables in the future by neutrinoless double beta decay and the radiative 
emission of neutrino pairs from excited atoms [24) . but we do not count them as observables in this paper. 
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Just from the degrees of freedom, one can fit all the 19 parameters using the full degrees of 
freedom in the model, in principle. However, the parameters we need to ht is hierarchical and 
some of the degrees are phases, and consequently, the hts do not necessarily done contrary to 
the simple number counting expectation. Actually, due to the hierarchy rriu/mt <C mg/m,- <C 
the complex freedom r2 should be almost consumed to £t the up-quark mass 
(roughly, the solution is obeyed by a cubic equation of r2, det 0). Thus, roughly 

speaking, to ht rriu, two degrees of freedom is consumed. As we have explained, the ht does 
not depend on R for small R very much, and then the phase of 5 is not an active freedom due 
to the rephasing covariance described above. Therefore, the active number of parameters to 
ht the quark masses and mixings are reduced. As a result, for the case of small R, down-quark 
mass rud is not fully ht as we will see later, though R ~ 0(1) {vr ~ 10^^ GeV) solution in 
type II can reproduce m^, as it was known in literature [IT] . 


4.2 Method of fit 


Before going to the detailed description of the ht process, we explain the equations to solve 
the quark masses which we can hnd quickly from the expression. 


rM„ = Me + = Ml/‘^ 

1 - r2 


3 -|- r2 A' — 6i?l -|- {K — 6Riy + D1 


1-^2 


18R - 26 


M]0. (4.12) 


As we have already noted, due to the hierarchy mu/rrit me/rriT-, the matrix in the bracket 
in the above equation has to be nearly singular. Therefore, 


1 — r2 Aj — QR Si \/(Aj — QR)^ -|- D 

3 + r2 ISR - 2(5 


(4.13) 


for one of i = 1,2,3, where Aj are eigenvalues of the matrix K as was given previously. By 
perturbation from the above relation, the up-quark mass rriu can be easily ht. One can also 
easily solve the equation by r: 


+ ml + ml) = Tr MuM^ 


Tr 



8 + r2 
1 - r2 




8+ r2 
1-^2 



(4.14) 


and thus, roughly speaking, the top quark mass mt can be easily ht by using r. As we have 
explained, we will try to ht the masses and mixings in quark sector by inputting the masses 
and mixings in lepton sector. The parameters to ht are now 


'fTT'di '^usi ^cbi ^ubi JkM- 

The function is dehned as 

.2 _ (Xi “ Xi)^ 

^ / j ''1 ’ 

^ cr; 

't * 


(4.15) 


(4.16) 
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Fixed Values 




rriu 

3.961 X 10-^ GeV 

m 


Parameters to ht 

mt 

71.0883 GeV 

izg 


rric 

0.1930 ± 0.025 GeV 

1251 

rrie 

3.585 X 10-4 GeV 

ra 


rrid 

(9.316 ±3.8) X 10-4 GeV 

m 

rrif. 

7.5639 X 10-2 GeV 

1251 


rris 

(1.76702 ± 0.5) X 10-2 GeV 

125 ] 

rrir 

1.3146 GeV 

izg 


mb 

(0.9898 ± 0.03) GeV 

125 ] 


7.54 X 10-® eV2 

0 


Vus 

0.224 ±0.002 

m 


2.4 X 10-3 eV^ 

0 


Vcb 

(3.7 ±0.13) X 10-2 

P! 

^12 

0.583 

0 


Vub 

(3.7 ±0.45) X 10-3 

PI 

6^23 

0.710, 7r/4 

0 


<5km 

1.18 ±0.2 

PI 

diz 

0.156 

0 





Table 1; The reference parameters used in the tit. 


where Xi and dj are the experimental measurements of the parameters and their standard 
deviations of errors. 

The values and uncertainties used in the ht are summarized in Table [U where we use the 
value of the quark mass in Table III of Ref. [25] (MSSM, tan/3 = 10) and calculate the quark 
mixing angle and KM phases at GUT scale from PDG data [26] by DR scheme. 

For the minimization, we use the Metropolis-Hastings (MH) algorithm [27] of the Markov 
Ghain Monte Garlo (MGMG) method. Since we expect that there exist many local minima, 
we perform the replica exchange MGMG (REMG) sampling [28] where several copies of the 
Markov chain with various temperature Tj are simultaneously simulated. Each chain performs 
MH scheme under a likelihood function dehned by exp[—x^/(2Tj)]. Every after performing 
some MG steps, the state of the chains are swapped by comparing the values of the likelihood. 
For ordinary MGMG, the state of the chain is often trapped at a local minimum and costs 
a large amount of time when we search a likelihood function with many local minima. For 
REMG, on the other hand, a replica with high temperature (> 1) has a flatter distribution 
function and the sampling point can easily move on to another local minimum. Thanks to the 
swapping algorithm, this transition can propagate to low-temperature replicas and therefore 
we can effectively sweep all the local minima. 

We should note, however, that it is still challenging to hnd the absolute global minimum 
of the likelihood function with many local deep minima. We have made much effort to hnd 
the minimum including changing the initial parameter set, tuning the replica temperatures as 
well as the width of the parameter jump in MGMG sampling. Nevertheless, there is always a 
possibility that the likelihood has another better minimum due to our limited hardware and 
calculation time. 

The hts are done assuming that Ml is negligible (type I: 5 = 0) or non-negligible (type II; 
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Figure 1; Left panel: R dependence of for type I of 023 = 0.710 (red dashed line), 7r/4 (blue 
dashed line) and type II of 023 = 0.710 (red solid line), 7r/4 (blue solid line). Right panel: The 
same figure as the left one but we transform R to the energy scale crvr. 


6 7^ 0)@ under two values of 023 (0.710, vr/4). In the MCMC, we fix the scale parameter R to 
precisely investigate the behavior against the scale parameter, and therefore, the function 
has seven (nine) free parameters for type I (type II) h| We adopt uniform prior for the phases 
(oe; *^2, ^3, <5pMNS) Arg((5)) from —TT to TT and the magnitude of the seesaw parameter 

|(5| < 500. Log flat prior is applied to mi over the range of 10 < — logio(mi/GeV) < 16. 

4.3 Fit results 

Fig. [1] shows R (and crVpi) dependence of values. For all cases, the global best £t is 
found at around vr ~ 0(10^^) GeV where the y^ value is 1.5-2 for type I and < 0.001 for 
type II. The y^ value of type I mainly stems from down quark mass (~ 0.44 MeV) which 
is smaller than the experimental value by ~ 1.2-1.3(T. For type I, another local minimum 
appears at vr ~ 0(10^®) GeV where large deviation is given by rrid (~ l.Sa smaller) and 
rris (~ 2cr larger). For type II, the y"^ curve becomes flat above vr ~ 0(10^'^) GeV and no 
prominent local minimum exists. As we expected in Section 3, vr 0(10^®) GeV is disfavored 
for type I while type II gives a y‘^ value of ~ 2-2.5 where the deviation is dominated by the 
down quark mass (~ 1.4cr smaller). Compared with the case of 023 = 0.710, although the 

® The parameter S is given in Ea. (l2.17p . and the classification of the y^ minimal solution is described in 
section 3. 

^ We also search the global best fit and the fit at i? = 0.001 allowing the neutrino oscillation parameters 
(012,013,023) and neutrino mass difference (ArUg^j, free within a range of 3a. We find that most of 

the parameters except for 023 converge around their central values and we decide to fix these reference values. 

For 023, on the other hand, although the 023 deviation does not significantly change the y^ values, its 
uncertainty is still large. Moreover, the 023 value has a correlation with dpMNS S't minimum as one can 

see in section 4.4. Therefore, in our fit, we adopt two values of 023 as the reference value, where one is the 
global fit of the three neutrino oscillation parameters [5] (023 = 0.71) and the other is the maximal mixing 
(023 = 7r/4). 
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difference is not obvions, 623 = 7r/4 gives slightly larger valnes and shifts the minima to 
lower energy scale. Tables |2] and [3] show the snmmary of the ht parameters and each 

global best ht {vr ~ 0(10^^) GeV), the local minimnm of type 1 {vr ~ 0(10^^) GeV) and 
R = 0.001 of type 11 {vr ~ 0(10^®) GeV). 

The nnmerical qnantities of the right-handed nentrino mass are also obtained from the ht. 
For type 11, we hnd that the sign combination of y/K (51,^2,53) becomes a permntation of 
(—1,1,1) at a low R region {R < 0.1). Thns, as investigated in section 3, the lightest right- 
handed nentrino mass is small compared with the other two eigenvalnes. For example, at R = 
0.001, the lightest mass is 6.9 (5.5) x 10® GeV while the other masses are 3.3 (3.5) x 10^^ GeV 
and 1.3 (1.2) x 10^® GeV for 623 = 0.71 (vr/d). In addition, the lightest right-handed nentrino 
mass is insensitive to R while the other two eigenvalnes of the nentrino mass behave as ~ 1/R 
in that small R region, which is also consistent with the analytical consideration in section 3. 
Aronnd the global best ht {vr ~ 10^® GeV), in contrast, the lightest nentrino mass grows and 
becomes nearly the same scale as the second largest nentrino mass. For example, at the global 
best ht R, the lightest nentrino mass is 2.4 (2.2) x 10^® GeV while other two eigenvalnes are 
4.6 (3.2) X lO^o GeV and 11.1 (6.9) x lO^^ GeV for 623 = 0.71 (vr/d). 

We assnme the normal hierarchy of the nentrino masses, rrii < m 2 < m^, in the above 
discussion. We also search for x^ minimum for the inverted hierarchy case, which can be done 
since the neutrino mass matrix is given as input in our formula. We hnd that the ht does not 
lead to a competitive result within the energy scale from 10^® GeV to 10^® GeV, which gives 

> 200 . 
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Type I 


Best fit 

Local minimum 

to 

CO 

0.710 

7 r/4 

0.710 

7 r/4 

R 

2.9853 

4.7315 

0.23713 

0.29006 


1.2200 

-0.15753 

0.97741 

1.5312 


-0.77358 

1.9305 

1.7110 

0.95362 

ar 

2.2100 

-2.1821 

-1.7127 

-2.3612 

0^2 

-2.8818 

2.8578 

-2.7599 

-2.8910 

Q!3 

-0.45528 

1.8924 

-1.1936 

-2.0470 

<5pMNS 

-0.80638 

1.9094 

-1.5819 

-2.3119 

logio(mi/GeV) 

-11.501 

-11.273 

-11.381 

-11.270 

rric (GeV) 

0.1955 

0.1999 

0.2153 

0.2098 

TUd (GeV) 

0.0004810 

0.0004270 

0.0004380 

0.0004416 

rris (GeV) 

0.01814 

0.01798 

0.02705 

0.02791 

mb (GeV) 

0.9916 

0.9885 

0.9789 

1.004 

Kc 

0.2242 

0.2241 

0.2243 

0.2245 


0.003703 

0.003674 

0.004025 

0.004011 

Vub 

0.03692 

0.03692 

0.03693 

0.03612 

Sku 

1.218 

1.189 

1.206 

1.190 

Pull 

rrtc 

0.099 

0.275 

0.893 

0.672 

rrid 

-1.186 

-1.328 

-1.299 

-1.290 

rris 

0.094 

0.061 

1.876 

2.048 

mb 

0.061 

-0.042 

-0.362 

0.464 

Vuc 

0.086 

0.033 

0.156 

0.229 

Vsb 

0.006 

-0.057 

0.722 

0.691 

Vub 

-0.065 

-0.064 

-0.051 

-0.675 

^KM 

0.190 

0.045 

0.128 

0.048 

r 

0.0138 

0.0144 

0.0279 

0.0285 

r2 

0.323 + 0.00618i 

0.311 - 0.0033B 

2.78 + 0.00190* 

2.77- 0.379* 

crVr (GeV) 

2.85 X 10^3 

1.76 X 10^3 

4.69 X 10^^ 

3.81 X 10^^ 


1.48 

1.85 

6.70 

7.51 


Table 2: The fit result for type I. Pull is defined hy (xi ~ Xi)l^i ^or each observable. 
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Type II 


Best fit 

R = 

0.001 

O23 

0.710 

7 r/4 

0.710 

7 r/4 

R 

5.1582 

7.2861 

0.001 

0.001 


-0.81968 

-1.8785 

-0.66648 

-0.25301 


-1.0015 

1.5169 

-2.8148 

2.8177 

Or 

1.4632 

2.8853 

-0.53961 

-0.84287 

O2 

-2.8481 

2.8306 

-2.8709 

-3.1146 

«3 

-0.30601 

1.1731 

-1.9809 

-2.8604 

(^PMNS 

-0.28115 

1.4476 

-2.3550 

-3.1131 

logio(mi/GeV) 

-11.592 

-11.467 

-11.207 

-11.173 

|5| 

75.056 

100.11 

15.545 

16.156 

Arg((5) 

0.11782 

0.11535 

0.43912 

0.51567 

rric (GeV) 

0.1931 

0.1929 

0.1978 

0.1989 

rrid (GeV) 

0.0009278 

0.0009309 

0.0004138 

0.0003936 

rris (GeV) 

0.01785 

0.01767 

0.01980 

0.02028 

rrih (GeV) 

0.9897 

0.9898 

0.9903 

0.9901 

Kc 

0.2240 

0.2240 

0.2240 

0.2241 

Vsb 

0.003698 

0.003698 

0.003765 

0.003724 

Vub 

0.03700 

0.03699 

0.03694 

0.03695 

Skm 

1.180 

1.180 

1.195 

1.160 

Pull 

rric 

0.004 

-0.005 

0.191 

0.236 

rrid 

-0.010 

-0.002 

-1.363 

-1.416 

nis 

0.036 

0.000 

0.426 

0.522 

mb 

-0.002 

0.000 

0.017 

0.010 

Vuc 

0.002 

-0.012 

-0.000 

0.027 

Vsb 

-0.005 

-0.005 

0.144 

0.054 

Vub 

-0.004 

-0.009 

-0.044 

-0.039 

^KM 

-0.001 

-0.001 

0.075 

-0.099 

r 

0.0140 

0.0138 

0.0230 

0.0233 

r2 

0.506 + 0.0252i 

0.502 + 0.00628* 

2.15 + 0.227* 

2.22 - 0.160* 

crVr (GeV) 

1.19 X 10^3 

0.861 X 10^3 

8.86 X 10^® 

9.22 X 10^® 


0.001 

0.0003 

2.10 

2.35 


Table 3: The fit result for type II. 
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^PMNS/^ 

Figure 2: The preferred region of 5pmns and 623 given by /g?. The best fit is shown by the red star. 
The red, blue, gray area represents la, 2a and 3a region respectively 


4.4 Prediction of PMNS phase in the neutrino oscillation 


The PMNS phase (5 pmns is one of the interesting predictions of our htting results. Though 
the accuracy is not still sufficient, the on-going experiments [3] already give a tendency of the 
parameter region. The long-baseline accelerator experiments and its global £t [29] suggest that 
623 larger than vr/d is slightly favored while the reactor experiments support small 6^23 (< 7r/4). 
Fig. 2 shows the global best fit by [5] which involves both the accelerator and the reactor 
experiment^. Although the wide range of the parameter region is still allowed, negative 
5pMNS with small ^23 (< 7r/4) is slightly favored. It is expected that the experiments will 
present more accurate ^pmns as well as 623 within a few years. In addition, accessibility of 
other future experiments is estimated and found to reach ten and several percent accuracy for 
6*23 and 10 to 40 percent accuracy for (5 pmns [SQ]- Therefore, it is interesting to compare the 
current experimental favored region with our £t prediction. 

For that purpose, we search for the global minimum for each hxed 6*23 and (5 pmnJ^ 
in the range of 0.675 < 6*23 < 0.835 and |5pmns| < allowing R free within a range of 


® The global best fit is recently updated by [Bj. The result does not significantly changed. 

® One might think that we can construct the same map by simply accumulating the MCMC samples by 
allowing 623 and i5pmns Fee. However, as mentioned in section 4.2, the function has a large number of deep 
local minimum and it is difficult to explore the detailed map under the full free parameter space. Hence, 
we decide to fix these two parameter for each MCMC simulation. 
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<5pMNs/^ 
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Figure 3: Left panel: The map of type I global best fit. The colored tiles show the magnitude 
of Middle panel: The map of type II global best fit. Right panel: The map of type II 
best fit at R = 0.001. We set the darkest blue panel to represent x^ > 25. 


0.0001 < R < 100 and also explore the same parameter region at i? = 0.001 {crUr ~ 10^® GeV) 
for Type II. Here we set the parameter range of 623 as 2a region of Ref. [5]. 

Fig. 3 shows the x^ map with respect to 5pmns and 623 for each case. The darker blue region 
represents higher x^ where the step size is given by the color bars beside the map. For type I 
global best £t (left panel), the figure shows that |(5pMNs/7r| < 0.3 is mildly disfavored when 623 
is large (> vr/d). At most of the points, the largest contribution to x^ comes from (smaller) 
rud with 1-2.5 c 7 deviation. The second largest contribution is from (larger) rris and (smaller) 
(5 km, which becomes conspicuous at the disfavored region and reach l-1.5a. Interestingly, 
the x^ asymmetry on the (5pmns sign is mainly derived from the (5km deviation, the center of 
which is located at around (5pmns/7I' ~ 0.1. If (5km is more precisely measured in the future, 
this asymmetry will be stressed, which leads a stronger prediction on the sign of the PMNS 
phase. 

For type II, the x^ map of the global best ht (middle panel) shows that x^ is small 
enough under any value of 6*23 and (5 pmns- In contrast to the type I case, larger x^ appears 
at |(5pMNs/7i'| > 0.5 in a large 623 (> 0.75) region, implying that the ht becomes an interesting 
prediction if 623 is large as suggested by the current accelerator experiments. The largest 
contribution to x^ is from rrid (0-1.2(T smaller) and the second largest contribution is from rris 
(0-0.7(T smaller) and the other contributions are negligible. We also check the x^ behavior at 
large 6*23 (0.835 < 6*23 < 0.88) under fixed (5pMNs/7r = —0.5 and find that the x^ value remains 
less than 3.0 for all the 623 region. We expect that more accurate measurement against the 
observables used in our ht will resolve the degeneracy of the x^ map in the near future. 

For type II of i? = 0.001 (right panel), though the shape of the x^ map is similar to the 
type I best ht, it seems more symmetric with (JpMNS-axis and its disfavored region is more 
obvious. The x^ value within |(5pMNs/7r| < 0.5 steeply increases as 6*23 becomes large and 
reaches ~ 80 at (5 pmns/ 7I' ~ 0. (In the figure, we set the darkest blue panel to represent 
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Figure 4: The relation between 5pmns Bind 623 under the two-zero texture assumption: = 

{M.u)i 2 = 0. The solid line is drawn using the center value of 613 : sin^is = 0.156, and the dashed 
lines are drawn to present the 2a range of 9 i 3 . For the mass squared difference and 612 , we use the 
center values of the global analysis 


> 25.) Particularly, ^pmns/tt ~ —0.5 with a large 623 (> 0.81) is strongly disfavored in this 
phenomenologically preferred case. Stronger upper bound on 623 can be found when |(5pmns| 
is small and reaches 6*23 < 0.73 at |5pMNs/7r| ~ 0. At a point with a small (< 5) in our map, 
the largest contribution to comes from (1 - 1.5cr smaller) and (0- 1.5cr larger). While 
the contribution does not strongly depend on 6*23, (5pmnS) fhe deviation becomes large 
at the disfavored region (|(5pMNs/7r| < 0.5 with a large 623 ) and reaches ~ 2.5a. Still, these 
contributions are subdominant at the disfavored region in which (larger) rric and (smaller) (5 km 
give 4-6(j and 3-4.5a deviations, respectively. In addition, (larger) mb and (larger) Vub also 
contribute to with 3-4a deviation around |(5 pmns/ 7’'| ~ 0 at a large 6*23 (> 7r/4). Similarly 
to the type I case, the (5km deviation is slightly (5pMNS-axis asymmetric, which is centered at 
|(5pMNs/7r| ~ 0.1. However, this asymmetry is canceled by the contribution from rUc, mb which 
are located at around |(5pMNs/7r| ~ —0.1. Hence, we cannot conclude which sign of ^pmns is 
preferred only from our £t. 

Even though it is not easy to explain how the prediction of the phase is obtained explicitly, 
one can understand the results qualitatively, especially in the case of i? = 0.001 for type H. 
Relating to the result that the £t of the down quark mass is smaller than the observation, 
(1,1) and (1,2) elements of are smaller than the other elements in the ht result. In fact. 
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in the special case with two-zero texture of the neutrino mass matrix, = 0 , 

one obtains a relation among the mixing angles, mass squared ratio, and PMNS phase EB. 
Under the two-zero texture assumption, we obtain [32l [33] 


COS(5pMNS — 


COS ^1^13 


V 

4 sin^ 6 >i 3 (^1 ^12 j sin 29i2 tan 6*23 


In Fig.4, we plot the relation between (5 pmns and ^23 in the assumption. Surely, those two 
elements are not exactly zero in the hts, and the relation provide a guide to understand the 
£t results for the prediction of the PMNS phase depending on the mixing angle. In fact, the 
points far from the two-zero texture, such as ((5 pmns5^23) ~ ( 0 , 0 . 8 ), are disfavored in our fits 
in the case of i? = 0.001 for type II. There is also such a tendency in the hts in type I, though 
the detail of the hts is more complicate. In the type II global best ht, on the other hand, the 
(1,1) and (1,2) elements are not small for the observed values to be tuned. Therefore, the 
values become slightly larger near the points from the two-zero texture. We remark that the 
smallness of ( 1 , 1 ) and ( 1 , 2 ) elements relates to the suppression of the proton decay amplitudes 
[34] . It is important to measure the PMNS phase and the 23-mixing angle accurately by the 
long-baseline neutrino oscillation experiments in order to survey the structure of the Yukawa 
matrices in the S'O(IO) model. 


5 Implication of the scale vr 

We have demonstrated the fitting of the fermion masses and mixings to make the properties 
of solution depending on vr conspicuously. As we have described, the key to distinguish the 
scale Vr in the model is the combination of the precise measurements of 23-mixing angles and 
PMNS phases. Surely, the hts are shown by making the analyses, and the prediction of the 
PMNS phase is modest by referring the values in the parameter space, though the hts will 
be more crucial when the parameters are measured more accurately in the future. Not only 
the numerical values of the observables but also the havor structure in the Yukawa matrices 
are diherent for each solution as we have displayed in the previous section. The detail of the 
havor structure can be applied to the low energy phenomena, which can depends on the scale 
Vr. In this section, we study the application, which can depends on the havor structure in the 
several cases of the model. 

5.1 Lepton flavor violation 

The htting of the fermion masses and mixings in the minimal S'O(IO) model hxes the matrices 
of the Yukawa couplings. The most interesting application of this feature is to calculate 
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the flavor violation in the SUSY model. In fact, the Dirac and Majorana nentrino mass 
matrices are determined (as a fnnction of vn scale), and the flavor violation is predictable. 
The importance of the scale vr can be also illnstrated in the scheme where the SUSY breaking 
mass matrices are diagonal as boundary conditions and the flavor violation is generated only 
from the Yukawa interaction via RGE running. As an application of the solution of the mass 
and mixing fitting, we perform the calculation of lepton flavor violation of p —?■ ey process 
by assuming the minimal flavor violation, which can be a major sensitive probe to And a 
footprint. 

There are two kinds of amplitudes of /i —)■ ey decay, Ar : fiR ^ ei'j and Ar : —)■ ej?y, 

and the branching fraction of /i —)■ ey is 

Br(/i ey) oc (5.1) 

Roughly speaking, Ar arises from the chargino loop contribution and the off-diagonal elements 
in the left-handed slepton mass matrix contribute to it. The right-handed amplitude arises 
from neutralino loop diagrams : (1) Bino-Higgsino diagram, (2) Bino-Bino diagram. The off- 
diagonal elements in the right-handed charged lepton mass matrix contribute to the diagram 
(1), and those of both left- and right-handed slepton mass matrices can contributes to the 
diagram (2). 

The sources of the flavor violation are the following: 

1. Dirac Yukawa coupling, Eu^Hu- 

The right-handed neutrino and the Higgs doublet propagate in the loop, and the 
off-diagonal elements of the left-handed slepton mass matrix are generated j35]. In the 
minimal S'O(IO) model, the Dirac neutrino Yukawa coupling and the heaviest right- 
handed neutrino mass are determined (irrespective of vr), and thus, this contribution is 
predictable. 

2. Majorana coupling, HAr. 

The SU{2)r triplet Ar and the left-handed lepton doublet i. propagate in the loop, and 
the off-diagonal elements of the left-handed slepton mass matrix are generated. In the 
triplet seesaw, when the SU{2)r triplet contribution is in the sub-eV range, the mass 
of the Ar is about 10^^ GeV, and this contribution will be important to the lepton 
flavor violation. Since the coupling matrix is unified to the 126 Higgs coupling, the 
Majorana coupling matrix is determined by the fermion mass and mixing fitting, up to 
the overall mixing parameter. The coupling matrix / is obtained as the original matrix 
f multiplied by the Higgs mixing as given in section 2, and thus, this contribution is 
not fully determined by the fit of the fermion masses and mixings. Since the matrix is 
obtained up to the over all factor, the ratio of the contributions to /i —)■ ey, r —)■ /ly and 
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r ^ 67 is predictable if the contributions from the other sources are small. Because the 
Higgs mixing cannot be larger than 1, the matrix elements of the Majorana coupling 
have lower bounds in the simplest model. However, if the 54 Higgs is adopted, the 
SU{2)l triplet from the 126 is mixed with the one in 54, and then, the SU{2)l triplet 
mixing is also multiplied, and the lower bound from this source is mo del-dependent. 

3. couplings. 

If the SU{2)ji symmetry or U{1)r x U{1)b-l symmetry is broken much below the GUT 
scale, the fields of A^“ and/or A^ can be light and they propagate in the loop and the 
off-diagonal elements of the right-handed charged-lepton mass matrix is generated. This 
contribution is related to the Higgs spectrum and model-dependent, and thus, it is not 
very predictable. However, if vr is ~ 0(10^^) GeV, this contribution is not small and 
need to be considered. Since it contributes to the FCNG from the right-handed sleptons, 
flavor violating decay is generated from the neutralino loop diagram. 

4. e^u'^Hc coupling. 

In the language of SU (5) GUT, 10 ■ 10 ■ includes this coupling. The right-handed up- 
type quarks and the colored Higgs fields Hq propagate in the loop, and the off-diagonal 
elements of the right-handed charged-lepton mass matrix are generated. In the other 
sources, the mass of the fields which propagate in the loop is supposed to be of the order 
of 10^^ GeV, while the colored Higgs mass is 0(10^®) GeV, and thus, this contribution 
is not very important, compared to the others in the current context of ^(^(lO) model. 

In general, those Yukawa coupling Yij to the charged-leptons can induce off-diagonal ele¬ 
ments in the slepton mass matrices by RGE in the form of 

E >■' (5-2) 

Here mg is a universal scalar mass, Aq is a universal scalar trilinear coupling, M* is a cutoff 
scale, Mx is the mass of a heavy field which propagate in the loop, and G is a group weight 
factor. In particular, the off-diagonal elements induced by the Dirac neutrino Yukawa coupling 
Yy can be roughly written as 

= - A + Al) In (5.3) 

where Yy is given in the basis where the right-handed neutrino mass matrix is diagonal, and 
is an eigenmass of k-i\i generation right-handed neutrino. It is important to note that the 
gluino and squark masses should be heavy due to the LHG results, and thus, the amount of 
the induced FCNG becomes less in the universal SUSY breaking models. The discovery of the 
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125 GeV Higgs boson also pushes up the squark masses. If the mass of the squarks and gluino 
are 0(10) TeV, it is hard to extract the off-diagonal elements from the flavor data. However, 
if the squark and gluino masses are about 2 TeV, the scalar trilinear coupling Aq has to be 
large (~ 5 TeV) to obtain the Higgs mass to be 125 GeV, and then, the off-diagonal elements 
are generated (even if mo is small) and the FGNGs are induced slightly. Therefore, the SUSY 
contribution can be consistent with the experimental results of many of the FGNG processes, 
but a slight excess can be observed in a process whose amplitude can have an enhancement 
factor. We remark that the circumstances are changed from the literatures a few years ago. 

Other main circumstance changed by the LHG experiments is the SUSY mass spectrum: 
The gluino and squark masses are bounded from below. The observation of Bg —)■ /i’*'//” to 
be consistent with the SM prediction also bounds the parameters [36]. Though the bound of 
the slepton masses, which is important for the lepton flavor violating decays, is not directly 
related to the gluino and squark mass bound, the slepton masses depend on the gluino and 
squark mass bounds if the universal SUSY breaking is assumed. The bounds depend on the 
SUSY breaking scenario. Since the anomaly mediation is an infrared phenomenon, the slepton 
spectrum is still less bounded compared to the squark masses, even in the GUT model£§. We 
calculate the lepton flavor violation using the Yukawa couplings by htting the fermion masses 
and mixings in the minimal S'O(IO) model assuming the universality to make the number of 
SUSY breaking parameters less, and thus, the slepton mass spectrum depends on the gluino 
mass bound indirectly. We choose the unihed gaugino mass mi /2 = 800 GeV to satisfy the 
current gluino bound and to be covered the next run at the LHG. We vary the universal scalar 
mass mo, which is important to the slepton masses. The scalar trilinear coupling Aq is chosen 
to make the Higgs mass to be 125 GeV, depending on mo. The ratio of the vev of up- and 
down-type Higgs bosons, tan /9, is an important parameter because the decay width of /i —>■ ey 
is roughly proportional to tan^ jS. We choose tan /3 = 10 in the calculation. 

We show the results of the numerical calculation in four cases given by the minimum 
hts, where the right-handed neutrino masses and the Dirac neutrino Yukawa couplings are 
given as follows in the basis where the charged lepton and the right-handed neutrino mass 
matrices are real/positive diagonal. Using unphysical phase freedom in the SM, we make the 
(f, 1) components to be real. 

1. type H, crVr = 8.86 x 10^® GeV (i? = 0.001) 

Mr^ = 6.9 X 10® GeV, Mr^ = 3.3 x 10^^ GeV, Mr^ = 1.2 x 10^® GeV, (5.4) 

In fact, if we take the known anomaly of the muon g — 2 seriously, the universal SUSY breaking is not very 
favored after the LHC results bound the gluino and squark masses, and the mixed modulus-anomaly mediation 
is favored in the unification scenario m- In that scenario, the branching fraction of /r —^ ey becomes larger 
and the flavor non-universality (as we will describe in Section 10) is needed to satisfy the current bound. 
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(5.5) 


/ 0.000111 0.000203+ 0.000217i 0.00888 + 0.00372i \ 
0.000440 0.0308 - 0.0248 i 0.0426 + 0.0013 i , 

\ 0.00607 -0.0276 - 0.0069 i 0.990 -0.278 i / 

2. type II, crVr = 1.19 x 10^^ GeV {R = 5.1582) (the best fit) 


Mr^ = 2.4 X 10^° GeV, 


Mr^ = 4.6 X 10^° GeV, Mr^ = 1.1 x 10^^ GeV, 


(5.6) 


Y, 


0.00301 -0.00147 - 0.00394 i -0.0031 + 0.0134 i \ 

0.0211 0.00023 -0.0284 i -0.0073 + 0.0359 i , 

0.103 -0.044 - 0.131 i -0.324 + 0.428 i / 


3. type I, crVr = 4.69 x 10^^ GeV {R = 0.23713) 


(5.7) 


Mr^ = 6.6 X 10® GeV, Mr^ = 1.3 x 10^^ GeV, 


Mr^ = 5.5 X 10^® GeV, 


(5.8) 




0.0000957 0.000056 - 0.000189 i 
0.000222 0.0279 - 0.0253 i 

0.00400 -0.0081 + 0.0177i 


0.00852 - 0.00459 i \ 
0.0418 - 0.0143* , 

-0.02-1.Hi / 


(5.9) 


4. type I, CrVr = 2.85 x 10^® GeV {R = 2.9853) (the best ht in type I) 


Mr^ = 6.8 X 10® 


H = 


GeV, Mr^ = 3.6 X 10^^ GeV, 

= 3.0 X 10^2 GeV, 

(5.10) 

0.000443 

0.00236 + 0.00016* 

0.0024 + 0.0109* \ 


0.00302 

-0.00648 - 0.00218* 

0.0105 + 0.0413* , 

(5.11) 

0.0208 

0.100 - 0.0004i 

0.004 + 0.527* / 



As explained, the effects from e^e^A~^~, couplings are absent in the case 1, and 

in type I solutions (3 and 4), the contribution from the Majorana coupling is absent. The 
contribution from the left-handed Majorana coupling has an ambiguity from the Higgs mixings, 
and here we assume the mixing to make the contributions to the flavor violation maximal. In 
the hgure, the behavior of the plot in case 3 looks different. This is because the (3,3) element 
of the Dirac neutrino Yukawa coupling is a bit larger than others, and the stop masses are 
reduced indirectly due to the RGE evolutions, and a larger value of the A-term is needed 
to reproduce 125 GeV Higgs mass compared to the others. (The (3,3) element in case 1 is 
also large, but the effect is not large since the 3rd right-handed neutrino is heavy and soon 
decoupled in the RGE evolution.) The induced off-diagonal elements in the slepton mass 
matrices are proportional to 3mQ -|- Aq, and thus, the ttiq dependence of Br(/i —)■ ey) looks 
different between the smaller and the larger rriQ. As noted, the value of Aq is chosen to make 
the Higgs mass 125 GeV, the detail numerical values depend on the Higgs mass. We note 
that due to a large value of Aq (~ 4-6 TeV), the stau becomes LSP for ttiq < 600 GeV. The 
branching fraction of the case 2 is the largest (in the region satisfying the current experimental 
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Figure 5: The branching ratio /x —)• 67 for type II at crvr = 8.86 x 10^® GeV (R = 0.001) (red 
solid line), type II at the global minimum (crvr = 1.19 x 10^^ GeV) (red dashed line), type I at 
the local minimum (crvr = 4.69 x 10^“^ GeVj (blue solid line), and type I at the global minimum 
(crvr = 2.85 X 10^^ GeVj (blue dashed line). We set 623 = 0.71. The gray dashed line shows the 
current experimental bound, Br(p ej) < 4.2 x 10“^^ j38jl . 


bound). This is because the (1,1) and (1,2) elements of the Dirac neutrino Yukawa coupling 
are larger than the other cases, relating to the result of that the proper size of the down quark 
mass £t in the case 2 . 

Relating to the vr dependence, we here comment on the electron EDM. We assume Aq 
and the Higgsino mass p to be real because otherwise the EDM becomes much larger than 
the current experimental bound. Even in the case, if there are e-r flavor changes (with CP 
phases) for both left- and right-handed slepton mass matrices, the Bino-Bino diagram can hit 
the At-V (1 — prUr tan (3 for the electron EDM amplitudes, and the value of the EDM can be more 
than the current experimental bound for slepton masses to be ~ 1 TeV. Surely, if there is only 
left-handed flavor change such as the case of crVr = 10^® GeV, the electron EDM is tiny as 
long as Aq and p are real. Therefore, if the scale of vr is much less than the unification scale, 
and the coupling can generate the right-handed slepton flavor changing, the electron 

EDM can be large in the S'O(IO) model, and the electron EDM can be the probe of the scale 
Vr. 

In the numerical calculation in the above cases, the electron EDM becomes large except 
for the case 1 for vr ~ 10^® GeV. In the case 2, the current bound of the electron EDM, 
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\de\ < 8.7 X 10“^® e • cm [39], can be satisfied for mo > 2.5 TeV, and then, we obtain Br(/i —>■ 
ej) < 1.1 X 10“^^. In the case 4 (the best fit in type I), we need mg > 2 TeV, and then 
Br(/i —>■ ey) < 2.3 x 10“^^. In the case 1, the electron EDM is tiny (< 10“^® e • cm) and 
mo ~ 1 TeV is still allowed, which can provide the boundary value of the current p —)■ ey 
bound. 

5.2 Proton decay 

As we have explained, the flavor structure is different in each solution, which influences the 
proton decay amplitudes depending on the scenarios of seesaw mechanism and the scale vji. 
The calculation of the proton decay amplitudes in various scenarios has been done in literatures 
[40] . and we do not repeat it in this paper. However, in the earliest works, they sometimes 
neglect right-handed proton decay operators to calculate the decay rate, and the importance 
of the flavor structure (especially for the (1,1) and (1,2) elements in the Yukawa matrices) for 
natural suppression of the proton decay amplitudes jSj is missing. We stress that the flavor 
structure is important to suppress the proton decay amplitudes in this model, and claim the 
importance of the large vji solution. We can show that the proton decay amplitudes can be 
suppressed near the line shown in Fig.4 in the solution. The detail will be discussed in a 
separate paper [41] . 

5.3 Discussion on the modification of Yukawa couplings by thresh¬ 
old corrections 

The quark masses and mixings can receive large corrections from the low energy threshold 
effects. This inversely means that the present procedure in the S'O(IO) model (inputting 
the lepton parameters and outputting quark masses and mixings) provides a prediction to 
holomorphic Yukawa couplings for quarks at the unification scale, and the deviations from 
the observed values of quark masses can be buried by the threshold effects [laialllll. This 
feature can be applied to suppress the dimension-five operators in the SUSY GUT models. 

As given in the previous section, in the case R 0(1) (vji ~ 0(10^^) GeV), the down 
quark mass (and all the other masses and mixings) can be fully fit. In the solution, the 
(1,1) element of the 126 Higgs coupling (in the basis where is diagonal) is as large as the 
down quark Yukawa coupling. On the other hand, the solutions for i? <C 1, the (1,1) element 
can be much less than the down quark Yukawa coupling, which has a merit to the natural 
suppression of the proton decay amplitudes though the fit of the down quark mass deviates 
from the observation at around 2cr. 

In the SUSY models, in addition to the corrections from the RGBs, the finite loop cor¬ 
rections arise due to the SUSY breaking. It is famous that the bottom quark mass can be 
corrected largely by the finite corrections of the non-holomorphic term, qd^H*. Actually, the 
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quark mass matrices, which are connected to the GUT scale via RGBs, can be modified by the 
loop correction via the SUSY breaking around TeV scale, and in general, all the quark masses 
and mixings can be different from them via the RGE evolution. However, the corrections 
from the loop corrections are severely restricted from the flavor physics data, such as K-K, 
B-B mixings, b —)■ sy, and Bg^ —)■ and there is not so much room to be available. For 

example, if the bottom quark mass is modified, the 6 —?■ sy and Bg^ —)■ processes can be 

also modified from the SM predictions. If the second generation masses are modified, the K-K 
and/or D-D mixings are modified. One may adjust them to the flavor data in general even 
if there are large corrections in all the quark masses and mixings, but it is not an attractive 
situation. Among the masses and mixings, only the first generation masses can be modified 
without a major contradiction to the flavor data, and thus, it is suitable to make the minimal 
modification to the fit results. 

The finite correction of the down-type quark masses from the gluino loop diagram (ne¬ 
glecting the effect from flavor mixings) can be written as 



where rrig is a gluino mass, F{x) = logx/(x — 1), and 9lr is a left-right mixing angle. 


sin 6] 


A'-jVd - jJ^mdi tan ^ 


(5.13) 


Thus, if the A term is proportional to the Yukawa coupling as in the minimal supergravity 
model {A^^Vd — Arridi), ^^e gluino contribution to the fermion mass is flavor universal : 
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(5.14) 


The chargino contribution can break the flavor universality to the mass corrections, but the 
corrections to down and strange quark are still universal. In order to correct the fit of the 
down quark masses and mixings, we need to break the universality of the corrections, and 
therefore, the non-proportional term of the trilinear coupling is needed. 

In the minimal supergravity model, the scalar potential contains the trilinear coupling 
originated from the mixed term between the visible and hidden sectors, W = HAis + fUhid, and 
thus, the scalar trilinear couplings are proportional to the Yukawa couplings. In the string- 
inspired models, the Yukawa coupling can depends on moduli, r, and the non-proportional 
term cau be mdueefl The Yukawa matrix and the derivative of the Yukawa matrix are not 

In general, denoting the cubic coupling of the superfield in the superpotential as YjjK^i^j^K, one 
obtains the scalar trilinear coupling as 


Aijk = Aq(Yijk cdrYuKir)). 


( 5 . 15 ) 










necessarily simultaneously diagonalized, and therefore, the non-proportional term can have 
complete different structure compared to the fermion mass matrix. As a result, the (1,1) 
element of the non-proportional part of the A term (in the basis of the Yukawa matrix is 
diagonal) is not necessarily hierarchically small, and it can easily modify the hrst generation 
fermion masses. If the gluino and squark masses are around 2 TeV, we estimate that cos P 
is sub-GeV to modify the down quark mass. Suppose that the all the components are roughly 
same order of the trilinear coupling, the corrections to the other elements of the mass matrix 
are also of the order of 1 MeV, which can be negligible for the generation mixings and the 
FCNC processes. 

If the down quark mass is modihed via the trilinear scalar coupling, the up quark mass 
can be also modihed. Therefore, the (1,1) element of the trilinear coupling for the up-type 
quarks should be small, which is possible within a freedom of the model parameter. Because 
there are only two Higgs coupling matrices in the minimal model, the (1,1) element of the 
trilinear coupling for the charged lepton remains in that case, and thus, the electron mass is 
also modihed. The correction to the electron mass is induced by the Bino component of the 
neutralino, and the modihcation is not large. Instead, there can have freedom to tune the loop 
correction of the up quark mass and cancel the tree-level one to make the physical up quark 
mass tiny. 

The scalar trilinear triplet Higgs couplings can induce dimension-four proton decay opera¬ 
tors via the triplet Higgs scalar exchange, and the (1,1) element of the coupling can make the 
excess of the proton decay amplitudes after double gaugino dressing. However, to generate the 
dimension-four operators, the holomorphic bilinear mass term of the triplet Higgs scalars is 
needed. Besides, the other elements of the trilinear scalar coupling matrix are free and can be 
small. Thus, the induced amplitudes can be suppressed within the freedom of the parameters. 

As the implication to the low energy observables, the bounds from the neutron EDM need 
to be considered, as in the case of the minimal flavor violation, and the non-proportional term 
of the trilinear coupling should be hermitian (in the basis where the quark mass matrix is 
real/positive diagonal). Because the coupling matrix is symmetric due to the 50(10) algebra, 
the phases of the components in the non-proportional term are severely restricted. The non¬ 
proportional term can be a new source of the flavor violation via the RGE evolution of the 
scalar masses, and the /i —)■ ey process will the most stringent constraint on it. As described 
before, the non-proportional term is not necessarily have a hierarchical structure parallel to 
the Yukawa hierarchy, and within the freedom, the contributions to /i ^ ey can be switched 
off. Rather, choosing the size of (1,2)/(2,1) elements of the trilinear coupling, the size of 
the /i —ey amplitude can be adjusted, and the production studied in the previous sections 
becomes just a reference value. However, the qualitative statement relating to the vr scale 
dependence is expected to be kept. 

We note that the charge/color breaking (GGB) vacua may appear for the large (1,1) element 
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of the scalar trilinear coupling, and the electroweak symmetry breaking vacua may become 
metastable [la He], similarly to the models in which the fermion masses and mixings are 
modihed by the loop corrections [H] . Surely, if the lifetime of our vacuum can be larger than 
the age of the universe m, the model is still valid. The bounds from the tunneling processes 
of the vacua should be analyzed model-independently, and the bounds are beyond the scope 
of this paper. We just mention that the vacua may be unstable, but long-lived as in the other 
analyses. 

6 Conclusions and discussions 

We study the minimal 5'O(10) model in which quarks and leptons couple to only 10 and 126 
Higgs representations. In the minimal Yukawa coupling, the up quark mass is not necessarily 
small, and as a consequence, the relation mem^rrir — mdrrtsmh is not obvious for a proper size 
of Cabibbo angle [H] . However, the quark and lepton masses can be fully htted (without such 
a clear reason) by using the full parameter freedom even in the case of the minimal Yukawa 
coupling. In fact, as we have described, in the case of R r'U 0(1) (i.e. vr ~ 10^^ GeV), the 
freedom is active to £t all of the fermion masses and mixings, which is consistent with the 
number counting of the parameters. 

Even though the fitting of fermion masses and mixings is just a choice of parameters, the 
minimal S'O(IO) is still attractive since the model is predictive. As was seen, the scale of vr 
is important for the minimum hts, and it is related to an intermediate scale determined by 
the S'O(IO) symmetry breaking vacua. Although we perform the search of the minimum 
by using the fermion masses at the GUT scale in the minimal SUSY boundary conditions, the 
qualitative properties of the solutions described in Section 3 are also applied to the case of 
non-SUSY boundary conditions. In non-SUSY model, the intermediate scale is rather favored 
for the gauge coupling unihcation, contrary to the SUSY models, and the constraints from the 
dimension-hve proton decays and the flavor physics are of course absent. To say it inversely, 
the predicted fermion mass matrices cannot be probed nor applied to the low energy physics 
in the non-SUSY model. In the sense to probe the GUT scale physics, the SUSY scenario of 
the predictive minimal Yukawa coupling in the S'O(IO) model is still attractive and worth to 
be chased. As it was explained, the solution for vr ~ 10^® GeV is interesting to apply even 
though the threshold corrections are needed to £t the center value of the down quark mass. 

As methodology, physicists often assume the minimality of a model, not only because the 
model is predictive due to the reduced number of the parameters, but also because the proto¬ 
typal model can contain the essential features to study the phenomenology of the objects. In 
general, it is obvious that the predictivity can be obtained by reducing the number of param¬ 
eters. It is important to dissect the property of the objects to see whether it only originates 
from the minimality or whether it comes from any other fundamental mechanism theoretically. 
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In reality, the minimality at the GUT scale physics is hard to be verihed experimentally, and 
one can always claim that the corrections can be indnced since the reduced Planck scale is 
just two digit above the GUT scale. Therefore, the concern one has to care about is whether 
the prediction is stable under the perturbation of the minimality assumption in GUT models, 
and we should investigate both the speciality of the prediction from the minimality assump¬ 
tion and the essential features which remain predictable even after the model is extended. As 
described, the major prediction of the minimal S'O(IO) model (in which the Yukawa coupling 
to the fermions is minimal) is that the Dirac and Majorana neutrino mass matrices are fully 
determined. 

Sometimes, the Higgs contents are also restricted to be minimal, i.e. 10 + 126 + 126 + 210, 
and the number of parameters in the Higgs superpotential is made to be minimal in the 
S'O(IO) model. In the SUSY version of the minimal Higgs contents, the fermion data httings 
require a typical S'O(IO) symmetry breaking vacuum, and some decomposed particles from 
the Higgs multiplets appear at the intermediate scale to obtain the proper size of neutrino 
masses. The gauge coupling evolution is severely affected by the presence of the particles at the 
intermediate scales, and the gauge unihcation suggests a non-SUSY (or split-SUSY) version of 
the model da ESI [IT]. In non-SUSY scenario, however, it is impossible to communicate with 
the predicted fermion mass matrices by using the low energy phenomena, such as the various 
modes of baryon number violating processes and lepton flavor violation, which are induced 
by the Yukawa interactions. We, thus, take a stance to probe the footprint of the GUT scale 
physics by using the SUSY scenario, in which the whole picture of the predicted mass matrices 
can be influenced to the low energy physics. In order to have the stance, we do not employ 
the minimality of the Higgs contents. As we have described frequently in this paper, we dehne 
the minimal S'0( 10) model as that based on the minimality of the Yukawa interaction, which 
is fully predictive to hud the structure of the fermion mass matrices. 

The important parameters in SUSY GUTs, which can affect low energy physics, are 126 
vev vr, SU{2)r triplet Higgs mass Ma, the lightest colored Higgs mass and the heavy 

gauge boson mass Mq. The parameters, vr, Ma, and Mq, are functions of the pa¬ 

rameters in the Higgs potential and particle spectrum (depending on the consistency of the 
gauge coupling evolutions) deductively, and one can discuss how the parameters can depend 
on S'O(IO) breaking vacua and the particle spectrum [19]. Instead of specifying the Higgs po¬ 
tential and deriving the scale vr from the certain number of parameters in the Higgs potential, 
we choose them to be free parameters as methodology because there are plenty of choices to 
extend the Higgs potential. In fact, vr is an important scale to describe the GUT symmetry 
breaking and the neutrino mass spectrum, and it is an open question how the right-handed 
Majorana neutrino is generated in GUTs and the size of the light neutrino masses are ob¬ 
tained. As we have seen, the £t of the fermion masses and mixings depends on vr. Actually, 
the naive scale of the right-handed neutrino masses to reproduce the observation is less than 
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10^^ GeV, which is hierarchically smaller than the GUT scale. In other words, the neutrino 
masses are a bit heavier if the GUT scale, 10^® GeV, is considered to be a fundamental scale, 
and explaining the hierarchy is one of the open questions in GUTs. In the solution of the 
minimal S'O(IO) model featured in this paper, the nearly singular Majorana mass matrix is 
contained for a large scale of vr instead of the smallness of the each component of the matrix. 
The nearly singular Majorana matrix is one of the solution to explain the neutrino masses 
keeping the symmetry breaking scale to be 10^® GeV, because the inverse of the Majorana 
mass matrix is included the type I seesaw formulae. In that sense, this solution can provide 
an approach to the fundamental question, and it is interesting to research the prediction of 
the model. 

In SU{5) GUT, the Dirac neutrino mass matrix is completely independent to the other 
charged fermions. In the unihed model scenario, the hierarchical forms of the Dirac and 
Majorana neutrino mass matrices are often assumed. The Dirac neutrino mass is generated 
by the couplings to the up-type Higgs, and in a certain model, the hierarchy of the Dirac 
neutrino mass is assumed to be similar to the up-type quarks. On the other hand, the predicted 
Dirac neutrino Yukawa coupling in the minimal S'O(IO) model is less hierarchical rather than 
the up-type quarks, and the hierarchy rather resembles to the down-type quark or charged- 
leptons. This is simply because the up-type quark mass hierarchy in this model is realized 
by a cancellation between two Yukawa matrices. Even in the non-minimal models, in which 
120 Higgs is added to the Yukawa couplings and the (1,1) elements of the Yukawa couplings 
are assumed to be small to suppress proton decay amplitudes in an S'O(IO) model, the Dirac 
neutrino Yukawa coupling also resembles to the down-type quarks and the charged-leptons 
and the size of the induced FGNGs is similar to the minimal model |33]. In that sense, the 
size of the induced FGNGs is predictable in the renormalizable S'0( 10) models. This feature 
provides an important implication to the process of /i —>■ ey in SUSY models, and it can be 
tested by the experiment. 

As described in the text, the electron EDM is one of the major implications of the vji 
dependence in the S'O(IO) model. If the scale vr is less than the unihcation scale, the flavor 
violation is generated in both left- and right-handed slepton mass matrices, and the electron 
EDM is enlarged. If the source of the flavor violation is hermitian, the phase of the EDM 
amplitude is canceled. In the minimal model, the Yukawa coupling is symmetric due to 
the S'O(IO) algebra, and the components are not real numbers to reproduce the KM phase. 
Therefore, the electron EDM can be induced in general for the lower scale vr. On the other 
hand, the Yukawa coupling to the 120 representation is anti-symmetric, and thus, there can 
be room to generate the hermitian flavor violation in the non-minimal models. Gonsequently, 
the electron EDM can be a major probe to distinguish the models. 

Among the quark-lepton phenomena directly induced by the Yukawa interaction, the GP 
violation in the neutrino oscillations is the last piece to be discovered. In fact, the ongoing 
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experiments have already given the hint that CP is violated in the oscillations, and it is 
expected that the CP phase (PMNS phase) is measured in the near future. Remarkably, 
the minimal S'O(IO) model has predicted the proper size of neutrino 13-mixing angle (before 
its measurement), which is consistent with the current experimental results. We show the 
prediction of the PMNS phase in the minimal S'O(IO) model, in which there are disfavored 
region depending on the 23-mixing. Now the prediction of the PMNS phase of the model 
is being challenged to the experimental results, including the precise measurement of the 
23-mixing. 


A Appendix : Square root matrix 


A square root matrix of A, denoted as \/A, is dehned to be matrices which satisfy 

In our numerical practical reason, we assume that the eigenvalues of a 3 x 3 matrix A, 
-^ 1 , 2 , 3 ; are not degenerate. Then, it can be diagonalized by a matrix V: 


Ai 


A = V 


As 




(A.l) 


A.? 


The square root matrix can be expressed as 


i-v/W 

VA = v{ ±y% 


v-\ 


±\/A^ 


(A. 2 ) 


The double-signs are in no particular order (INPO), and there are 8 -fold matrices for the 
square root of a 3 x 3 matrix. 

Dehning matrices Aj as 


0 


0 


Ai = V 

we obtain 


0 


A2 = V 


0 


1 I As = R 

0 


0 


A — AiAi -|- A 2 A 2 -l- A 3 A 3 . 

One can easily obtain the following relations: 

Ai -|- A2 -|- A3 = 1, AjAj = Aidij. 

By solving the simultaneous equation 

1 1 1 \ / A, \ /I 

Ai A2 A3 1 j A2 1 — I A 

A? Ai Ai / V A3 / V 


R-\ (A.3) 

(A.4) 

(A.5) 

(A. 6 ) 


33 


we obtain 


Al = 

(a.-a.)(a,-A 3 ) P (V + AjM + VA 3 I), 

(A.7) 

A 2 = 

(A 2 -A,)(A 3 -As)P (Ai + AsM + AiAjI), 

(A. 8 ) 

A 3 = 

(A3-A0(A3-V)P (A. + VM + A.VI), 

(A.9) 


Using the matrices Aj, the sqnare root matrix can be written as 

\/A = ^Siv^Ai, (A.IO) 

i 

where s, stands for the sign of each root of eigenvalne, Sj = ± 1 . 
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